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Abstract. Using partially twisted boundary conditions we compute the K — > 7r semi-leptonic form factors 



in the range of momentum transfers < q < q n 



{mK—fn-n) in lattice QCD with Nf = 2+1 dynamical 



flavours. In this way we are able to determine / + T (0) without any interpolation in the momentum transfer, 
thus eliminating one source of systematic error. This study confirms our earlier phenomenological ansatz 
for the strange quark mass dependence of the scalar form factor. We identify and estimate potentially 
significant NNLO effects in the chiral expansion that guides the extrapolation of the data to the physical 
point. Our main result is f+ n {0) = 0.9599(34)(1^)(14), where the first error is statistical, the second error 
is due to the uncertainties in the chiral extrapolation of the lattice data and the last error is an estimate 
of potential discretisation effects. 

PACS. 11. 15. Ha Lattice gauge theory - 11.30.Rd Chiral symmetries - 12.15.Ff Quark and lepton masses 
and mixing - 12.38.Gc Lattice QCD calculations - 12.39.Fe Chiral Lagrangians 



1 Introduction 



The last five years have seen tremendous progress in high 
precision calculations from first principles of the K — > tt 
semi-leptonic vector form factor at vanishing momentum 
transfer, f+*(0). The results are combined with precise 
experimental measurements of \f+*(0)V us to obtain the 
CKM matrix element V us . It has been shown in many nu- 
merical studies [IHZ] that the observable f+*(0) can be 
computed with sub-percent level of precision on the lat- 
tice. All common efforts take advantage of the fact that the 
lattice prediction is necessarily exactly unity in the SU(3) 
limit in which m 2 = m 2 K . Thus leading lattice errors on 
the form factor are proportional to the mass difference 
and apply to the difference of the form factor from unity 
f+ w (0) — 1 = Af + $2 , where is the analytically known 
0(m^^m 2 K ,m 2 ) ) contribution to the form factor in 5(7(3) 
chiral perturbation theory ■ The precision of current 
lattice computations is well represented by the result of 



the RBC-UKQCD collaboration [6], 

Af = -0.0129(33) stat .(34) q2>x (14) a 



(1) 



The first error is statistical, the second is due in approx- 
imately equal parts to the chiral extrapolation and to an 
interpolation in the momentum transfer q 2 and the final 
error is due to cut-off effects. 

The primary purpose of this work is to remove com- 
pletely the systematic error in /^(O) due to the ^-inter- 
polation at the lightest quark masses. In addition we study 
in more detail the ansatz for the chiral extrapolation which 
was used to obtain the above result ([1]). To this end we 
build on our previous work on the g 2 -dependence of me- 
son form factors and stochastic quark propagators [TTJHLl] . 
The chiral extrapolation error above includes a contribu- 
tion from using a simulated strange quark which is a little 
heavier than the physical one. Here we confirm this esti- 
mate with new simulation results which allow us to inter- 
polate directly in the valence strange quark mass. We also 
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discuss phenomenological fits based on NLO SU(2) [T5] 
and on SU(3) [3 [9] chiral perturbation theory for the 
vector form factor and we discuss potentially significant 
NNLO effects. 

Our results reinforce the RBC-UKQCD collaboration's 
computation of f+*(0) which constitutes a crucial input 
to precision tests of the Standard Model and for con- 
straining its possible extensions [TJ] . The results presented 
should be interpreted as an intermediate step towards 
the extension of the RBC-UKQCD collaboration's effort 
towards simulations at smaller quark masses and at a 
smaller lattice spacing which are under way. 

After a brief discussion of the computational setup and 
a comparison of our current strategy with the conventional 
one, we summarise and discuss our results. 

2 Strategy 

The matrix element of the vector current between initial 
and final states consisting of pseudo-scalar mesons Pi and 
Pf respectively, is in general decomposed into two invari- 
ant form factors which parameterise non-perturbative ef- 
fects, 

(Pf(pf)K\p t ( Pt )) = 

f+ Pf (q 2 )(Pi +Pf) fl + f P ' Pf (q 2 )(p,-Pf), , ^ 

where q = pf — pi is the momentum transfer. For K — > n 
scmi-lcptonic decays is the weak current sj^u, Pi = K 
and Pf = 7r. We also introduce the scalar form factor 

/<f V) = f? (? 2 ) + / 2 /5V) , (3) 
m K - 

which satisfies f^ v (0) — f+*(0). In this paper we are pri- 
marily interested in computing f+™(q 2 ) for q 2 = but 
we also present new results for other values of q 2 . In con- 
trast to previous calculations of the form factor in lat- 
tice QCD [IHZ] we simulate directly at q 2 = by using 
partially twisted boundary conditions [T5J[IS]. One com- 
bines gauge field configurations generated with sea quarks 
obeying periodic boundary conditions with valence quarks 
with twisted boundary conditions [15H24] , i.e. the valence 
quarks satisfy 

ip(xk + L) = e ie "^(xk), (k = 1, 2,3) , (4) 

where ip is either a strange quark s or it represents one of 
the degenerate up or down quarks q. The dispersion rela- 
tion for a meson in a finite volume projected onto Fourier 
momentum pft then takes the form [THUD] , 

E = yj m 2 + (p FT + AO) 2 , (5) 

where m is the meson mass and where AO is the difference 
between the twist angles of the two valence quarks. 

In [TUJ we showed how this technique can be used to ex- 
tract meson transition matrix elements at arbitrary values 
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Fig. 1. Quark flow diagram for a 3pt function with initial and 
final states Pi and Pf, respectively. 

of the momentum transfer. In particular, for the matrix 
element with the valence quark flow diagram as in figure 
[I] with the initial and the final meson carrying 3- momenta 
Pi = Pft.i + Oi/L and pf = Pftj + f /L, respectively, 
the momentum transfer between the initial and the final 
state meson is 

q 2 = & -p f ) 2 = (Efa) - Ef(pf)) 2 ( Pl -p f ) 2 . (6) 

In our study it will be sufficient to set the twist of the spec- 
tator (anti-)quark q% to zero so that it satisfies periodic 
boundary conditions. 

As will become clear below, the matrix element in @ 
can be extracted from the time-dependence of combina- 
tions of Euclidean two- and three-point correlation func- 
tions which can be computed straightforwardly in lattice 
QCD. The two-point function is defined by 

C i (t,p i ) = T, x e i P^(O i (t,x)Oi(0,0)) 

, ,2 (7) 

= ^-(e-^ + e-^-*)), 

where i = 7r or K, and where Oi are pseudo-scalar interpo- 
lating operators for the corresponding mesons = q"/§q 
and Ok = sj^q and we assume that t and T — t (where T 
is the temporal extent of the lattice) are sufficiently large 
that the correlation function is dominated by the lightest 
state (i.e. the pion or kaon). The constants Zi are given 
by Zi = (Pi | Oj(0, 0) | ) . The three-point functions are 
defined by 

Cp!/p f (ti,t,tf,pi,pf) = 

= Z v J2j pH * , - x) e ipv *{OAt f ,x f )V tl (t,x)O}(t i ,0)) 

Xf,X 

= Zv ^E ( P f(Pf)\V„(0)\Pi(Pi)) 

x[e(t f -t)e- E ^-^- E ^f-^ 

+ CM 0(t- t/ ) e -^( T +'«-')- E /(*-'/)} , (8) 

where Pij is a pion or a kaon, is the vector current 
with flavour quantum numbers to allow the Pi — > Pf 
transition and we have defined Zf = ( | O/(0, 0)| Pf ). 
We have introduced the constant c M which is c.q = —I 
(time-direction) and Cj = +1 for i — 1,2,3. Again we as- 
sume that all the time intervals are sufficiently large for 
the lightest hadrons to give the dominant contribution. 



RBC and UKQCD Collaborations: K — > n form factors with reduced model dependence 



3 



The vector current renormalisation factor Zy can be 
obtained as follows. For illustration we take <t <tf < 
T/2, in which case Zy is defined by 



ci^'^ (ti,t, tf ,0,0) 



(9) 



In the numerator we use the function CV (t,p) — (t,p) — 
e~ B "( T_t ) where Z v and E n are determined from fits 
to Ctt (t, 0) and using eqn. (0 (for tf = T/2 it is natural in- 
stead to use Cir(t,p) = 7}C„(t,p) in ©). The superscript 
B in the denominator indicates that we take the bare (un- 
renormalised) current in the three-point function. 

In the following we drop the labels f$ and tf (since 
they are fixed) and we combine the two- and three- point 
functions into the ratios 



R 



[% iPf (Pi,Pf)=M< 



c Pi {tf,Pi)c Pf {tf,p,) ■■ 



Kr C l P f P f (t,Pi,p f ) 

N — '- 1 



/ C Pi (tf -t, Pi ) C Pf (t,p f ) C Pf (tf ,pf) 
Cp f (tf,Pf) V C Pf (tf-t,p f )C Pi (t, Pt )Cp t (t f , Pt ) ' 

(10) 

where N = AZy^jEiEf and the ratios are constructed 
such that 



/* P ' (q 2 )( Pl + p f ), + (q 2 )( Pl - p } ), , 



111) 



for a = 1, 3. For the ratios we use the naming convention 
of [TU] but we haven't made use of ratio 

Once these ratios have been computed for some choices 
for pi and pf while keeping q 2 constant (of course we 
are particularly interested in q 2 = 0) the form factors 
f+ v {q 2 ) and f^^(q 2 ) can be obtained as the solutions of 
the corresponding system of linear equations. 



3 Simulation parameters 

We simulate with Nf = 2+1 dynamical flavours generated 
with the Iwasaki gauge action [35J[2S] at (3 = 2.13, which 
corresponds to an inverse lattice spacing a -1 = 1.73 GcV 
(a = 0.114(2)fm) [HIUS], and the domain wall fermion 
action [291130] with a residual mass of am rcs = 0.00315(2) 
[2"T] . The simulated strange quark mass, am s = 0.04, is 
close to its physical value [55], and here we choose the 
RBC-UKQCD configurations with the bare light sea quark 
mass am q = 0.005, which corresponds to a pion mass of 
330MeV [2"71l2"5] . The calculations are performed on a lat- 
tice volume of 24 3 x 64 sites with the fifth dimension hav- 
ing an extension of 16 lattice points. More details can be 
found in [2"g] . 

We distinguish three sets of correlation functions as 
summarised in table [T] The correlation functions in set 

1 We did not generate data for Cp M p(t, p, p) \ q 2 =0 for P = 
ix,K from which the forward matrix elements (P|V^|P) rele- 
vant for the construction of R2 can be extracted. 





P4 


Z4PSs4 


Z4PSs3 




0.04 


0.04 


0.03 


# meas 


700 


1180 


1180 


# sources 


4 


8 


8 



Table 1. Some basic simulation parameters. 



P4 (point source, 4 positions of the source) are identical 
to those on which the RBC-UKQCD prediction for /f ""(O) 
at am q = 0.005 in [B] is based. The naming convention for 
the remaining two data sets, ZAPSsA and ZAPSs3, is mo- 
tivated as follows: spin-diluted Z(2) x Z(2) pH[12] noise 
source and point sink with strange quark mass am s = 0.04 
and am s = 0.03, respectively. The data set ZAPSsA with 
light quark mass am q = 0.005 corresponds to a unitary 
simulation point, i.e. where the sea and valence quark 
masses are the same, while set Z4PS's3 corresponds to 
a partially quenched parameter choice. For Z4PS*s3 and 
ZAPSsA we started the measurement chains for eight dif- 
ferent source time-slices. In each case we measured on 
every 40th trajectory and averaged the correlation func- 
tions over the chains into bins of eight. The correlation 
functions obtained in this way were computed with zero 
Fourier momentum and the momenta of the initial and/or 
final kaon/pion were induced by twisting one of the kaon's 
and/or pion's valence quarks. For each measurement we 
applied the full twist along one of the spatial directions. 
We changed this direction frequently as the measurements 
proceeded in order to reduce the correlations. Our choices 
for the twisting angles for ZAPSsA and ZAPSsi and the 
corresponding values of q 2 are summarised in table [5] In 
order to obtain q 2 =0 we make the following two choices 
of the twisting angles [TU] : 



'K 



2m „ 



and K\=LJ(^) 2 



i 2 and 6k = 



(12) 



As input to these formulae we have used the estimates 
for the central values of the kaon masses arriK = 0.2990 
(ZAPS3) and am K = 0.3328 {ZAPSA) and for the pion 
mass am^ = 0.1907 (for both datasets). These values had 
been determined from a previous study of the gauge field 
ensemble considered here. The momenta of the mesons are 
given by px = Ok/L and p^ = O^/L. In addition to the 
values of 9^ and 9k in eqn. (jT2"]) , propagators were gener- 
ated for other values of the twisting angle. In particular, 
for the kinematical situation where the kaon is at rest and 
the pion is moving due to the additional ad hoc twist an- 
gle 9 n = 1.600 we determined the corresponding values 
for 9k which yield the same q 2 also when the pion is at 
rest. The contractions of these propagators into two- and 
three-point functions allow a computation of f+ n (q 2 ) for 
additional values of the momentum transfer in the range 
from about q 2 = to g^ax ( c f- table |2J. 
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6 n 


0k 


g 2 /GeV 2 








am s = 0.04 




2.682 


4.68 



0.0002(2) 
0.0004(3) 


0.9758(44) 


0.9758(44) 


-0.0997(93) 


2.129 





0.0216(2) 


0.9898(34) 


0.9975(42) 


-0.081(17) 


1.600 



0.0381(2) 
2.792 0.0382(2) 


1.0030(20) 


1.0213(32) 


-0.108(11) 








0.0607(2) 


1.0185(15) 






am s = 0.03 


2.682 





-0.0192(3) 


0.9677(49) 


0.9613(41) 


-0.054(14) 




2.129 


3.337 0.0001(5) 
0.0001(3) 


0.9867(30) 


0.9867(30) 


-0.0771(77) 


1.600 



0.0149(3) 
2.509 0.0150(4) 


0.9986(21) 


1.0066(27) 


-0.0852(96) 








0.0352(4) 


1.0124(10) 







Table 2. Table of twisting angles used in this study, together 
with the corresponding values of q 2 and the results for the form 
factors. 



4 Data analysis 

The results for the two- and three-point correlation func- 
tions for P4, Z4PSs3 and ZAPSsA were analysed using 
the jack-knife as well as the boot-strap procedure. In all 
cases covariance matrices for the correlation functions and 
ratios R\ and R3 were generated in order to be used for the 
fits (frozen covariance matrix). Unfreezing the covariance 
matrix, i.e. using the covariance matrix computed individ- 
ually for each jack-knife or boot-strap sample destabilised 
the fits. We interpret this as a reflection of the fact that 
we have an insufficient set of measurements and that the 
fluctuations of the covariance matrix are therefore large. 
We found that the results we get with the frozen covari- 
ance matrix agree within (similar) errors with the results 
we would get when neglecting any correlations. 

We found that for the spatial component of the vector 
current in the three-point function, i.e. for and R3, 
the plateaus are of comparable quality - in the analysis we 
opted to use R^. For the time-component however in the 
cases where only one of the initial and final states carries 
a twist the quality of the ratios varies strongly. Here we 
decided to use R { ° ] for the case where only the pion carries 

the twist and R$ in all other cases. 

Table [2] summarises the kinematical points which we 
analysed. The kaon masses for the full statistics of ZAPSs3 
and ZAPSA turn out to be arriK = 0.2987(4) and amx = 
0.3327(4), respectively and the pion mass in both cases 
is anin = 0.1903(4^3- From the table we see that there 
are degeneracies in q 2 , i.e. we have data for the same q 2 
but from three-point functions with different kinematical 



While agreeing within errors, the central values differ 
slightly from those quoted in [28] because the number of mea- 
surements and the measurement techniques differ. 



parameters for the kaon and pion. In the analysis we first 
compute the corresponding ratios R\ and R3 and then 
average them as indicated by the degeneracy. 



5 Extrapolation models 

In a previous analysis [5J we used results for the form fac- 
tors at the lattice Fourier modes for four pion masses and 
one value of the (unitary) strange quark mass which post 
facto turned out to be slightly heavier than the physical 
one. At each simulated light quark mass we first deter- 
mined two estimates for /i f "'(0) = /^"(O), one from an 
interpolation in q 2 with a polc-ansatz, 

/+"*' (0)| pole 



fo^il 2 ) I pole 



l-q 2 /M 2 



(13) 



and one from an interpolation of fo(q 2 ) with a 2nd order 
polynomial, /f 7r (0)| 2n d (cf. table IV in [6]). We based the 
estimate of the systematic error due to the phenomeno- 
logical interpolation on the difference between the two re- 
sults. 

The final central value was determined from a global fit 
ansatz incorporating pole dominance, the NLO expression 
for /f^O) in [8], 



/f 7r (0)|NLO = l + / 2 (/,m^m^) 



(14) 



(/ is the pion decay constant), and modelling higher order 
contributions. The ansatz is 



ml) 2 (A +A 1 (m 2 K + mD) 



q 2 /(M + M 1 (m 2 K 



where we use the N 



ml)) 

(15) 

/ = l -1- 1 expression for f 2 , par- 
tially quenched in the strange quark, as determined in [5] . 
Since the Kaon mass appears explicitly, after fitting (|15ll 
to our lattice data, any values for m n and can be 
inserted into eqn. (15) to obtain a value for / ( f 7r (g 2 ). 
Hence, by inserting the physical values for m n and mx, 
the strange quark mass is automatically corrected to its 
physical value. The data was well described with 



A Q = -0.34(9)GeV" 
M = 0.94(10)GeV, 



Ai = 0.28(12)GeW 6 
Mi = 0.54(18)GeV _1 



(16) 



In the meantime an expansion for f+*(0) in SU (2) chi- 
ral perturbation theory has been derived Q3] . In particular 
in view of the slow convergence of SU(3) chiral perturba- 
tion theory observed for some quantities (cf. e.g. [28]) it 
seems useful to compare the present extrapolation strat- 
egy to one incorporating this new formula. Similarly to 
the case of SU(3) chiral perturbation theory we use the 
ansatz, 



F+(l 



c 2 ml 



c 4 ?r4) 



l-g 2 /(M + Mim; 



(17) 
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Fig. 2. Summary of simulation results of /<) (0). The black 
circles and the (dashed ) pole interpolation correspond to the 
results of [S] while the results represented by the left- and right 
pointing arrows, respectively, correspond to the results of this 
work for am s = 0.04 and am 3 = 0.03. The red and blue curve 
represent the result from the global fit in [B], once for m^ 03 



where L = jg^jp log {j^TJ with the normalisation scale 
\x and where in comparison to the original work we have 
added an additional term proportional to m%. Note that 
the parameters in this fit ansatz depend on the value of 
the strange quark mass. 

Reordering the chiral expansion: We do not really 
know how light the quarks must be for the chiral expansion 
at a fixed order to represent the mass dependence of phys- 
ical quantities to a given level of precision. In principle, 
(unrealistically expensive) lattice computations at very 
light quark masses, perhaps even lighter than the phys- 
ical ones, could answer this question. Our present calcu- 
lations however, involve masses in a regime where NNLO 
terms of the SU (3) expansion are non-negligible, and yet 
we have insufficient data to determine these terms and all 
the corresponding low-energy constants (LECs). As dis- 
cussed above, we therefore have to model the higher order 
contributions using an ansatz such as that in eqn. (1151) . 

We observe that the interchange symmetry {q 2 ) = 
f+ K (q 2 ) is held in the 5/7(3) expansion order by order 
and also in non-perturbative data to all orders. The SU(3) 
chiral expansion in terms of the unknown, but in principle 
unambiguous, LEC fo has this symmetry manifest in each 
term. However, we have the freedom to repartition terms 
of this expansion between different orders: for example to 
use an alternative expansion parameters / differing from 
fo beyond leading order 

f^(0) = l + f 2 (f,m 2 K ,ml) + ... . 

In fact, the NLO term f 2 is usually quoted as 
hifnt m K> m 7r) — —0.023, with the physical pion decay 
constant used in place of the unknown LEC fo- For SU(3) 
chiral perturbation theory to correspond to QCD when all 



terms are summed, this must be a passive transformation 
on the series making simultaneous but cancelling changes 
to different orders. There are several possible interpreta- 
tions of the f.„ which appears in f 2 {f ni m|, m 2 ,) that are 
consistent with QCD in the asymptotic light mass region. 

First, we may consider /„■ entering the expression to be 
simply an estimate of the true LEC fo, and a rigorous error 
analysis must allow for to vary over a range of values. 
In order that the NNLO form is unchanged, and that the 
interchange symmetry remains manifest in our ansatz, we 
choose this first interpretation and vary the value of fo 
entering our expressions to estimate a systematic error. 

Second, a mass dependent expression for 
f7r{foi m ud,nT's) may enter f 2 and the correction term 

£nnlo = h{fo, m 2 K, ml) - f 2 (fir(fo, m ud , m s ), m 2 K , ml) , 

introduced at NNLO, passively repartitioning the series 
and possibly improving convergence. We note that this 
difference term now breaks the symmetry under n <H> K 
interchange, and the form of NNLO terms is changed to 
compensate the alternate expansion parameter. 

A third but less helpful interpretation would be to con- 
sider adjusting the series to expand in terms of the physi- 
cal decay constant f n ~ 131 MeV. Passive transformation 
requires that a correction term at NNLO be introduced 
with the same dependence on simulated masses 

<5nnlo = f2(fo,m%,m%) - f 2 (U = 131MeV, m 2 K , ml). 

This simply relabels the term involving the free parame- 
ter fo as NNLO. If fo is subsequently expanded in terms 
of fn and m^ t ph,y S , this amounts to using a phenomeno- 
logical estimate for fo — with the associated systematic 
uncertainty. 

Using f 2 (f = 131MeV, m 2 K , ml) but failing to adjust 
the forms appearing at NNLO is inconsistent with QCD 
in the true chiral limit as it actively changes the series. 
This discussion impacts our previous analysis in which we 
extrapolated the defect Af, and did not admit variation 
in fo 7^ 131 MeV, nor modified the form of our global 
fit ansatz at NNLO to admit breaking of the mass in- 
terchange symmetry. While we are concerned here with 
the context of the K — >• it decay, the considerations are 
applicable to practices used in the lattice study of many 
quantities if a chiral expansion parameters differing from 
fo is used to improve convergence of the truncated series. 

6 Results 

Cost of the simulation: In table [3] we compare the cost 
of the simulations for the two approaches to the compu- 
tation of f+*(0). For P4, for each quark mass am = 
am q , am s one normal and one extended propagator (cf. 
definition in [TO]), one twist (periodic boundary condi- 
tions), four positions of the point source on 175 configura- 
tions and 12 spin-colour inversions were necessary. Since 
with point sources a Fourier transformation between the 
source and the point of the current insertion can be per- 
formed almost for free, one can directly interpolate to q 2 = 
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am n props n g ?i src n conflg s-c N tot 

PA 2 x 2 x 1 x 4 x 175 x 12 = 33600 

ZAPSsA 2x 2 x 2 x 8 x 147 x 4 = 37632 



at no additional cost. This is not the case when using the 
noise source technique as for ZAPSsA and ZAPSsi. We 
achieved a similar precision for f+ v (0) at approximately 
the same total cost with the ^4PS'-source type, where for 
each propagator of mass am = am q or am = am s four 
spin-colour inversions are necessary [TT] for each choice of 
the twist angle. We note however that in general the qual- 
ity of plateaus is significantly enhanced when using the 
stochastic volume source technique. 

<7 2 -dependence of the form factor: The data gener- 
ated for this paper complements our previous data PA by 
a number of new points for f^^iq 2 ) in the range < q 2 ^ 
<7^ lax for two strange quark masses am s = 0.04 (unitary) 
and am s = 0.03 (partially quenched). The results are il- 
lustrated in the plot in figurc[5]by the red/blue right/left- 
pointing arrows, respectively. The new data points for 
am s = 0.04 nicely agree with both the pole dominance 
and polynomial fits (cf. eqn. (fi~3"|) in [B]) as can be seen in 
the following comparison: 

results for am q = 0.005, am s = 0.04 

/£*(<)) | P cie = 0.9774(35) 0, 

polynomial — 

0.9749(59) 0, 
/f 7r (0)|thiswork = 0.9757(44). 

In [BJ we used the spread /+ ""(O) | po i e - f+ w (0) Ipoiynomiai ~ 
0.0024 as an estimate of the systematic due to the phe- 
nomcnological ^-interpolation. As simulations move closer 
to the physical pion mass, the value of q 2 ^^ = (tuk — w w ) 2 
increases. Therefore the interpolation to q 2 = 0, which 
crucially depends on the high precision which one is able 
to achieve for the form factor at q 2 ^^, will be increasingly 
sensitive to the ansatz one makes. One therefore expects 
the systematic error due to the interpolation to increase. 
We emphasise that the approach advocated here entirely 
removes this uncertainty. 

Quark mass dependence: Inserting the unitary and 
partially quenched kaon mass which we simulated here 
together with the parameters in (|16[) into the phenomeno- 
logical ansatz (|l~5|) we can predict the form factor that 
is to be expected for am s = 0.03 and am s = 0.04 with 
am q = 0.005 as illustrated in terms of the blue (dot- 
dashed) and red (dashed) curve in figured Both curves are 
nicely compatible with the new blue and red data points, 
thus confirming that the ansatz parameterises the depen- 
dence of the form factor on a partially quenched strange 
quark well. 

Combining the data sets Pi, ZAPSsi and ZAPSsA 
and carrying out the global fit (| 15[) we update the previ- 
ous result f^{0) = 0.9644(33) — » f£*(0) = 0.9630(34) 
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Fig. 3. Result of global fit based on SU(3) chiral perturba- 
tion theory using /2(115MeV, mpc , m-n). The plot also shows 
a comparison with recent results for the form factor in chiral 
perturbation theory. 



(statistical errors only) at the physical point. The result 
of the global fit is also illustrated in figure [3] by the solid 
black line. 

The chiral extrapolation of the lattice data is well 
constrained by the natural hinge-point (0)\ mi< =vi n = 
1. As can be seen in figure [3J our data as well as the 
global SU(3) fit-ansatz (fl"5j) nicely approach this point 
for m^ —> mx- In contrast, in SU(2) chiral perturba- 
tion theory one expands the form factor around vanishing 
pion mass at a fixed strange quark mass [13) (in fact, all 
strange quark mass dependence resides in the low energy 
constants). The limit f^* '(0)| ms = m = 1 is not naturally 
implemented in this expansion. Given our experience with 
SU(2) fits to other pion and kaon observables in [25] such 
an expansion describes the lattice data reliably only be- 
low ps 400MeV. In contrast to our study in [5S] here 
we only have data for two values of the pion mass be- 
low this cut-off and extrapolations are therefore not well 
constrained. Alternatively one can include data at heav- 
ier pion masses. However, fits of acceptable quality can 
then only be obtained after adding an extra term (oc C4) 
to the expression in [T3]. Given these considerations, at 
this stage we refrain from presenting fit results based on 
SU(2) chiral perturbation theory. 

Estimates of systematic errors: The new data pre- 
sented here confirms the ansatz for the g 2 -interpolation 
for the smallest mass used in rcf. [BJ, i.e. am q — 0.005. 
Since q 2 lax increases as m q decreases, it is at this mass 
that 9 2 lax is the largest (and therefore furthest away from 
q 2 = 0) and hence the interpolation is the least con- 
strained. We are therefore confident that the pole ansatz 
previously used in fits to our data [B] describes the form 
factor data well also for all the other simulation param- 
eters where q^ax ^ s closer to the origin. The systematic 
error due to the interpolation can be safely removed from 
our final result. 
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Table 3. Cost comparison. 
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Fig. 4. Illustration of the dependence of the fit result (with 
the ansatz in eqn. (|15[) 1 on the choice of the decay constant. 
The horizontal red lines indicate our estimate of the resulting 
systematic uncertainty. 

As discussed in section [51 a potential source of sys- 
tematic error which to our knowledge has not been taken 
into account systematically in any previous computation 
of /+' ir (0) is the choice of the decay constant entering in 
the SU(3) NLO prediction for the form factor. Lacking a 
precise value of the decay constant in the chiral limit we 
repeated the global fit for the three choices / = 100, 115, 
131 Me V and found for the central values of the form- 
factor /f w (0) = 0.9556, 0.9599, 0.9630, respectively. In 
each case the fit was of very good quality. This is quite a 
sizeable variation in the central value which is illustrated 
in figure 0] We found that the choice of decay constant 
particularly changes the slope of /+' ir (0) with respect to 
in the region of small pion masses where we do not 
have data. In order to study the behaviour at NNLO more 
systematically, one can use the FORTRAN computer code 
written by Bijnens^l that can provide SU(3) NNLO terms 
in terms of numerical integration routines (see also |32j for 
form factor fits based on this code). Our experience from 
using the code is that for our limited set of lattice data 
points there are too many free parameters (low energy con- 
stants from the 0(p 4 ) and the 0(p 6 ) chiral Lagrangian) 
to be able to carry out reliable fits. Lacking a better an- 
alytical understanding of NNLO effects we prefer as the 
central value the one corresponding to / = 115 MeV. As 
an estimate for the uncertainty in the chiral extrapola- 
tion we use the interval defined by the result for /+' ir (0) 
as obtained when using / = lOOMeV and / = 131 MeV, 
respectively. Our result is therefore, 

jf-(0) = 0.9599(34)(i^)(14), (18) 

where the first line the first error is statistical, the sec- 
ond is due to the chiral extrapolation and the third error 
is an estimate of lattice cut-off effects for which we stay 
with the previous estimate in [5] . We note that the quoted 
uncertainty due to the chiral extrapolation covers the cen- 
tral value of the result which one obtains when extrapo- 
lating instead with ^^^{jn u ^m s \rr? K ^m^) in (15), i.e. 
with the decay constant as input that corresponds to each 
of our simulations points (c.f. [28] )• Adding all errors in 

3 the programs for f+(q 2 ) and f-(q 2 ) used in Ref. [31] are 
available on request from Johan Bijnens 



quadrature we obtain 

/**(0) = 0.960(+*). (19) 

We believe that the systematic error due to the chiral 
extrapolation discussed above is conservative but still a 
more rigid statement would be desirable. To this end a 
better understanding of the NNLO terms in the chiral 
expansion and additional simulation points at smaller pion 
masses are mandatory. 

7 Summary and Outlook 

In this paper we present new data for the kaon semi- 
leptonic decay form factors computed in lattice QCD with 
Nf = 2 + 1 flavours of dynamical fcrmions. Using (par- 
tially) twisted boundary conditions we were able to di- 
rectly simulate at the phenomenologically relevant kine- 
matical point q 2 = 0, at a very small additional computa- 
tional cost while at the same time making the computa- 
tion independent of any phenomenological ansatz for the 
interpolation in the momentum transfer. In this way one 
significant systematic uncertainty in the computation of 
the form factors has been successfully removed. 
We have reconsidered the estimates of the systematic un- 
certainties of our previous calculation presented in [B] and 
we present numerical evidence that our ansatz for the 
strange quark mass dependence of the form factor is un- 
der control by providing data at an additional partially 
quenched simulation point for the strange quark mass. 
Currently, chiral extrapolations of lattice results for the 
kaon scmi-lcptonic form factor are based on NLO chi- 
ral perturbation theory. We show that ambiguities in the 
parametrisation of the NLO expression can lead to addi- 
tional systematic effects which we include into our revised 
estimate of the systematic uncertainties. This ambiguity 
also applies to any other lattice computation of /^(O). 
We want to stress that the interpretation of lattice data 
for the K — > tt form factors would profit from the avail- 
ability of their expressions at NNLO in chiral perturbation 
theory in a more transparent form. 

The RBC-UKQCD collaboration is currently extending 
the set of data points presented here by simulations at 
lighter pion masses and at the same time at a finer lattice 
spacing. A combined analysis of all data is the next step in 
RBC-UKQCD's program of a precise computation of the 
K — > 7T form factors in Nf =2 + 1 flavour lattice QCD. 
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